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We present a versatile scheme for creating topological Bogoliubov excitations in weakly interacting 
bosonic systems. Our proposal relies on a background stationary field that consists of a Kagome 
vortex lattice, which breaks time-reversal symmetry and induces a periodic potential for Bogoliubov 
excitations. In analogy to the Haldane model, no external magnetic field or net flux is required. 
We construct a generic model based on the two-dimensional (2D) nonlinear Schrodinger equation 
and demonstrate the emergence of topological gaps crossed by chiral Bogoliubov edge modes. Our 
scheme can be realized in a wide variety of physical systems ranging from nonlinear optical systems 
to exciton-polariton condensates. 
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Introduction .—The quantum Hall effect is one of the 
most celebrated results of modern condensed matter 
physics [1]. The robustness of the Hall conductance can 
be traced back to the non-trivial topology of the underly¬ 
ing electronic band structure [2], which ensures the exis¬ 
tence of chiral edge modes and thus eliminates backscat- 
tering. Recently there was a surge of interest in the pos¬ 
sibility to exploit such topology to create chiral bosonic 
modes in driven-dissipative systems — with possible ap¬ 
plications for one-way transport of photons po- 

laritons HMS] , excitons nais], magnons nanH], and 
phonons [iniisn]. A common thread through these seem¬ 
ingly diverse ideas has been to induce topology by ex¬ 
ternal manipulations of a single-particle band structure, 
with interactions playing a negligible role. Exceptions 
from this non-interacting paradigm are proposals that 
combine strong interactions with externally induced ar¬ 
tificial gauge fields to create non-equilibrium analogs of 
bosonic fractional quantum Hall states [2T1124] . 

Here we take a new perspective and consider (bosonic) 
Bogoliubov excitations (“Bogoliubons”) where weak in¬ 
teractions induce a non-trivial topology [25] . We demon¬ 
strate that topological Bogoliubons naturally arise on 
top of a condensate that exhibits a lattice of vortex- 
antivortex pairs, with no net flux required. Interactions 
are key to harness the time-reversal symmetry breaking 
provided by the condensate vortices. Erom the point of 
view of Bogoliubov excitations, they generate non-trivial 
“hopping” phases which lead to an analog of the Haldane 
lattice model [26|. The periodic potential defining the 
corresponding lattice can be introduced as an external 
potential or generated by interactions with the underly¬ 
ing condensate. 

Although our scheme can be applied to any system de¬ 
scribed by a 2D nonlinear Schrodinger equation (Gross- 
Pitaevskii equation, or analog thereof), we focus on sys¬ 
tems of weakly interacting bosons that have a light com¬ 
ponent, where the required vortex lattice can readily be 
obtained from the interference of several coherent (op¬ 
tical) fields (see Eig. [^. In this setting, the phase¬ 



FIG. I. Top: Optical pumping scheme allowing to create a 
suitable condensate for topological Bogoliubons, with incident 
field composed of six equal-frequency plane-wave components 
with wave vectors kn as depicted by arrows and phases chosen 
as = 0 except for (f 4 = 27r/3 and (fe = —27tI3. Bottom: 
Intensity (left) and phase (right) pattern of the resulting field. 
Intensity maxima form a Kagome structure (white lines) with 
vortex-antivortex pairs located such that each hexagonal pla- 
quette is threaded by a flux 24>o = 47r, with smaller triangular 
plaquettes threaded by — 4>o. Although fluxes cancel out over 
the whole system, a well-defined chirality emerges, defined by 
the sign of 4>o. Bogoliubov excitations of the condensate ex¬ 
perience non-trivial fluxes due to interactions, leading to TR 
symmetry breaking and ultimately to topological states. 


imprinting mechanism through which non-trivial topol¬ 
ogy is ultimately granted is analogous to that proposed 
a few years ago in the context of optomechanical sys¬ 
tems [27j. The same mechanism was recently applied 
to create topological phonons using photons that remain 
trivial m- Here we demonstrate that vortex lattices can 
be exploited in a much broader variety of systems where, 
in contrast to previous works, topology does not emerge 
from the coupling to a different bosonic species, but is 
intrinsically enabled by interactions. 
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Our analysis starts with the 2D nonlinear Schrodinger 
equation commonly used to describe weakly nonlinear 
dispersive physical systems: 

= (wd(-iV) + Q!|V’(x,i)p + F(x)) V’(x,i), 

( 1 ) 

where '0(x, t) is a complex field (or “wavefunction”) de¬ 
scribing the amplitude and phase of a coherent field, and 
describes the dispersion of the system. Nonlin¬ 
earities are characterized by a parameter a which may be 
positive (repulsive/self-defocusing interactions) or nega¬ 
tive (attractive/self-focusing interactions). We also allow 
for a spatially dependent potential V (x) although, as we 
will show, this is only essential when a > 0. Equation 0 
is applicable to a wide variety of physical systems rang¬ 
ing from light propagation through Kerr nonlinear me¬ 
dia [28] to the dynamics of exciton-polariton systems [29] 
and Bose-Einstein condensates m- Below we outline a 
generic mechanism for creating topological Bogoliubov 
excitations, assuming that the underlying bosonic fields 
have a light component. We then discuss specific imple¬ 
mentations. 

Theoretical scheme .—The first and essential ingredi¬ 
ent of our scheme consists of a stationary field '0o(x, t) 
which exhibits vortex-antivortex pairs and intensity max¬ 
ima that form a Kagome lattice pattern, as illustrated in 
Fig-lH In bosonic systems with a photonic component, 
such a Kagome “vortex lattice” can be directly imprinted 
onto the field '0(x, t) using an optical coherent field (or 
“pump”) composed of six plane-wave components as de¬ 
picted in Eig.j^(see also Supplementary Material EH). 
Denoting by ^/^(x) the (bosonic) field operator associ¬ 
ated with the light component, such coherent pumping 
can be described by a Hamiltonian term of the form 
/ dx/(x)e-*“o*4'^ (x) + h.c., where /(x) and cjq respec¬ 
tively denote the pump spatial profile and frequency. In 
the mean-field limit where Eq. 0 applies, this results in 
an additional pumping term proportional to /(x)e“*“^°^ 
(see, e.g.. Ref. [32|). The required pumping field is read¬ 
ily obtainable with a spatial light modulator [33] or by 
passing light through a mask [34] . 

Remarkably, the vortex-antivortex pairs arising in the 
above Kagome vortex lattice are located in such a way 
that each elementary hexagonal plaquette is threaded by 
a flux 24>o = 47r (i.e., contains 2 vortices), while smaller 
triangular plaquettes are threaded by — 4>o (i.e., contain 
a single antivortex) [see Eig.[^. Since the incident pump¬ 
ing field has no net orbital angular momentum to trans¬ 
fer to the internal field, vortices and antivortices always 
come in pairs and fluxes cancel out over the whole sys¬ 
tem. Nevertheless, the stationary field '0o(x, t) breaks 
time-reversal symmetry, as revealed by the “chirality” 
defined by the sign of 4>o (which is set by the phases of 
the pumping field [31] and does not change under rota¬ 
tions or translations in the plane that leave the Kagome 
lattice of intensities invariant [35]). 


The fact that '0o(x, t) exhibits a staggered flux pat¬ 
tern with no net flux but a well-defined chirality is very 
reminiscent of Haldane’s seminal proposal for quantum 
Hall physics with no external magnetic field [26]. Here, 
however, the phase that an excitation would accumulate 
around any plaquette of the Kagome lattice (or around 
any loop in the plane) seems to be an integer multiple of 
27r (equivalent to no phase at all), which would naively in¬ 
dicate that TR symmetry is preserved and that topolog¬ 
ical states are out of reach. Remarkably, we demonstrate 
below that Bogoliubov excitations, through interactions^ 
do experience the TR symmetry breaking of the pump¬ 
ing field. This will allow us to access a non-equilibrium 
interaction-induced analog of Haldane’s model. 

To derive the spectrum of Bogoliubov excitations on 
top of '0o(x,t), we perform a linearization of Eq. . 0 - 
We define the slowly-varying field 0(x, t) = '0(x, 
where ujq is the pump frequency, and consider weak per¬ 
turbations (or “fluctuations”) of the form [36] 

(/)(x,t) = 0o(x) +I4(x)e“*“'^ -h u*(x)e*“'*^ (2) 


where 0o(x) denotes the rotating-frame counterpart of 
the stationary field '0o(x, t) = 0o(x)e“'^“^o^, u{yi) and u(x) 
are (in general complex) functions determining the spa¬ 
tial form of the fluctuations, and uj is the frequency of the 
perturbations which is kept complex in order to capture 
potential instabilities [36]. Plugging the above expres¬ 
sion into Eq. Q and neglecting second-order terms in 
u(x) and v(x) yields the Bogoliubov equation 


/ u>'(x) 

V -acpKx-f 


a</>o(x)^ 



where a;'(x) = u!d{—i'V) — Wq + 2q;|^o(x)|^ + V{x). This 
shows that Bogoliubov excitations experience both the 
intensity and phase pattern of the underlying field (x). 

In the case of repulsive interactions {a > 0), the ef¬ 
fective potential term 2a|0o(x)p tends to localize the 
excitations at the vortex/antivortex points of the lattice 
(see Eig.[^. To allow for topological states, however, it is 
crucial that excitations remain localized away from vor¬ 
tices/antivortices so as to pick up a phase when hopping 
around them. This can be achieved, e.g., by introducing 
an external potential V{x) with minima that coincide 
with the maxima of |0o(x)p. In the case of attractive in¬ 
teractions (a < 0), the effective potential 2(a|0o(x)p au¬ 
tomatically localizes the excitations at the desired points, 
thus obviating the need for any additional potential. 

According to Eqs. 0 and (|^, fluctuations described 
by u{x) and v{x) can be viewed as “particle-hole” analogs 
of each other: u{x) as a “particle” excitation with disper¬ 
sion co’'(x), and v{x) as a “hole” excitation with opposite 
dispersion — ci;'(x). The pump frequency ujq sets the rel¬ 
ative energy between the two. Interestingly, the minus 
sign present on the off-diagonal of Eq. 0 (a hallmark 
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of bosonic Bogoliubov excitations m) makes the corre¬ 
sponding matrix non-Hermitian. The stability of the sys¬ 
tem is then assessed by the imaginary part of uj: when 
Jm(co’) > 0 , fluctuations grow exponentially and 0 o(x) is 
an unstable solution of Eq. 0 - _ 

The off-diagonal coupling of Eq. 0 is crucial for our 
scheme: it harnesses the TR symmetry breaking encoded 
in the condensate phase and thus ultimately determines 
the size of the attainable topological gaps. The periodic 
potential 2 (a| 0 o(x)P + E(x) also plays an important role: 
(i) it generates band crossings in the dispersion of u (and 
similarly for v) at which gap openings are most easily 
obtained; and (ii) it ensures that Bogoliubov excitations 
“live” on the Kagome lattice pattern defined by the max¬ 
ima of | 0 o(x)p, which allows them to experience in an 
optimal way the fluxes induced by 0 o(x). 

Tight-binding analysis. —Intuitive understanding of 
the Bogoliubov excitation spectrum can be gained from 
the tight-binding limit of Eq. (|^, which follows from 
the low-energy behavior obtained in a deep potential 
2(a|0o(x)p + — cjo whose minima form a Kagome 

lattice. In the sector u (v), we denote the correspond¬ 
ing on-site energy by and assume that excitations 

“hop” with amplitude t(—t) between nearest-neighboring 
sites. The off-diagonal terms in Eq. 0 lead to couplings 
iLge^^^^ between sectors u and v on each site j. To 
model the phase pattern of (^o(x)^, we choose the phases 
ipj so that each triangular unit cell of the Kagome lattice 
carries a single vortex (see Eig. [^). 

The mechanism giving rise to TR symmetry breaking 
(and, in turn, to topological states) can be revealed most 
clearly in the limit 2Q ^ g where particles and holes {u 
and v) are coupled off-resonantly and weakly. In that 
case virtual transitions between particles and holes lead 
to a renormalization of the hopping amplitude and phase. 
In a perturbative picture [37], a particle-like Bogoliubov 
excitation can be turned into a hole-like excitation, hop 
to a neighboring site, and transform back into a particle- 
like excitation (see Eig. 0 )- Since the particle-hole cou¬ 
pling involves a site-dependent phase, excitations pick up 
an overall phase (“+” for particles and “—” for 

holes) [38], which leads to an effective hopping 



FIG. 2. (a) Bogoliubov spectrum in the off-resonant tight- 
binding limit (in strip geometry): the particle-like (positive) 
part of the spectrum exhibits three bands with topologi¬ 
cal gaps crossed by pairs of counter-propagating chiral edge 
states, reminiscent of a Kagome lattice model with staggered 
fluxes m- The red (blue) coloring indicates the degree of lo¬ 
calization at the lower (upper) edge. Parameters were chosen 
as Q/t = 6 and g/Q, = 2/3, in a strip geometry with peri¬ 
odic boundary conditions in the x-direction (see Supplemen¬ 
tal Material m for details), (b) Unit cell used in the tight- 
binding model, with phases for the on-site coupling chosen 
as (0, 27r/3, —27r/3) in the basis {A, B, C). (c) TR symmetry¬ 
breaking mechanism: off-resonant couplings between particle- 
and hole-like excitations (u and v) lead to virtual transitions 
which renormalize the hopping term t. 


the positive or negative part of the Bogoliubov spectrum 
would exhibit three bands touching at high-symmetry 
points of the Brillouin zone (Dirac cones at K and K’, 
and a quadratic band touching at T). The non-trivial 
flux ^ induced by interactions gaps out these degenera¬ 
cies, leading to topological bands with a non-zero Chern 
number [39] (see Eig. [^. To second order in the size 
of the topological gaps is given by m 


Ak = 


4(^7 — “t" t/2) 


Ar 


2{n - t){n + 2t)' 


(5) 


in counter-clockwise direction around the triangular unit 
cell. The corresponding flux per triangular plaquette is 
given by = ±3 arctan(\/ 35 ^/( 8 f 4 ^ — 5 ^)) 7 ^ 0. 

Equation 0 illustrates one of the key aspects of our 
proposal: the fact that interactions generically lead to a 
non-trivial flux $ (not a multiple of tt), which breaks TR 
symmetry from the point of view of Bogoliubov excita¬ 
tions. In the off-resonant tight-binding regime considered 
here, particles and holes are both individually described 
by an effective Kagome lattice model with staggered flux 
pattern as depicted in Eig.[^ If ^ was trivial (for ^ = 0), 


The Bogoliubov spectrum considered here in the off- 
resonant regime is stable, with only real eigenvalues. 

Exciton-polariton systems. —An example of realization 
of our scheme can be practically considered with exciton- 
polariton condensates in semiconductor microcavities, 
which are renowned for their strong nonlinearities [29] . 
The Bogoliubov spectrum of polariton condensates has 
been measured in photoluminescence [40] and four-wave 
mixing [41] experiments. The required amplitude and 
phase pattern of 0 o(x) can be imprinted optically as de¬ 
scribed above [42] . 
















4 


The fact that polariton-polariton interactions are re¬ 
pulsive implies that a (Kagome) periodic potential V (x) 
is required for our scheme (as discussed below Eq. 0). 
Such potential can be realized, e.g., by depositing thin 
metal films on the structure surface [43j [44] (as illus¬ 
trated in Fig. ^), by applying combinations of surface 
acoustic waves [45| or by etching micropillar arrays [46] . 
Alternatively, one can engineer the desired potential opti¬ 
cally m by taking advantage of the fact that polaritons 
have two possible circular polarizations (spin projections 
along the structure growth axis), and that polaritons 
with opposite spins typically interact attractively [481(49] . 
Specifically, one can consider Bogoliubov excitations with 
spin a = ±1 and use a component with spin —a to induce 
the desired potential. In practice, this can be achieved 
by pumping both components simultaneously with an el- 
liptically polarized incident field. 

Here we follow this all-optical approach and consider 
Bogoliubov excitations with spin a on top of a polariton 
condensate with spin components 0o,cr(x) and 0o,-cr(x) 
(defined as in Eq. 0). The relevant Bogoliubov equa¬ 
tion is then given by Eq. ^ with a = ai and E(x) = 
<^ 2 | 0 o,-cr(x)p, where ai and a 2 denote the strength of in¬ 
teractions between polaritons with parallel and opposite 
spins, respectively [31]. To provide a reliable estimate of 
the size of the topological gaps achievable in practice, we 
compute the stationary fields 0o,±cr(x) and the spectrum 
of Bogoliubov excitations with spin a in an exact and 
self-consistent way, without relying on any tight-binding 
approximation (see Supplemental Material [31]). Our re¬ 
sults are illustrated in Fig. 

In accordance with the above tight-binding analysis, 
the particle-like (positive) part of the low-energy Bogoli¬ 
ubov spectrum exhibits three Kagome-like bands with a 
clear topological gap between the lowest two bands m- 
With conservative parameters from existing experiments 
(taken from Ref. [51]), the topological gap reaches about 
0.06meV, which exceeds typical polariton linewidths of 
the order of jaeV [52]. Indeed, an important feature of 
exciton-polaritons is that they also exhibit dissipation, as 
light can escape through imperfect mirrors. Dissipation 
can be treated by introducing a decay term into the (par¬ 
ticle/hole) dispersion, ±co’c^(—iV) ^ ±cc;£^(—iV) —'i/(2r), 
where r is the polariton lifetime [36] . 

Discussion. —A variety of other bosonic systems would 
be suitable for the realization of our proposal. The non¬ 
linear Schrodinger equation 0 provides, in particular, 
a direct representation of Maxwell’s equations for light 
propagating through a nonlinear medium with polariza¬ 
tion transverse to the propagation direction ^ (where 
t ^ z \n Eq. 0 and X = (x,^) defines the lateral co¬ 
ordinates) [28]. In that case the potential E(x) can be 
induced optically. In particular, in nonlinear media with 
a strong electro-optic anisotropy, polarizations can be 
chosen such that a suitably polarized incident field ex¬ 
periences negligible self-nonlinearity and induces a sta- 
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FIG. 3. (a) Bogoliubov spectrum of exciton-polariton systems 
(in strip geometry), exhibiting similar features as the tight- 
binding spectrum of Fig.[^, except for the upper topological 
gap which vanishes due to the modified dispersion. Although 
particle- and hole-like excitations overlap in energy (around 
Ct; = 0), all eigenvalues are real, indicating that the system is 
stable. Parameters were taken from Ref. m, giving a low- 
energy topological gap of about 0.06meV [31]. (b) Practical 
realization in semiconductor microcavities composed of dis¬ 
tributed Bragg reflectors (DBRs) [43], with periodic poten¬ 
tial E(x) induced, e.g., via metal surface patterning (other 
mechanisms are discussed in the text), (c) Realization of the 
effective potential E(x) in a nonlinear optical system using 
coupled waveguides ESI. 


ble potential E(x) for optical fields with opposite (or¬ 
thogonal) polarization, which experience a strong non¬ 
linearity [54H56] . Alternatively, E(x) can be realized 
by refractive-index modulation, e.g., in arrays of coupled 
waveguides [57] [58] as depicted in Fig. [^. In this con¬ 
text, nonlinearities > 0.2cm“^ with significantly 

weaker losses < 0.04cm“^ have been demonstrated using 
fused silica waveguides m- 

Equation 0 is also relevant for the description of cold 
atom condensates, providing an alternative to existing 
schemes for topological states. Optical flux lattices (in 
which optically dressed atoms both experience fluxes and 
a periodic potential) have been proposed and realized in 
this context uni EB, with a variety of other available 
methods for optical phase imprinting [62H65] . The stir¬ 
ring of a cold atom condensate by rotating laser beams is 
also well known for realizing Abrikosov-type lattices [66] . 
Although this would also be suitable for generating topo¬ 
logical Bogoliubov modes (via Eq. (|^), the transfer of 
total angular momentum by stirring would imprint a net 
global flux, unlike in our scheme. In that case the Bo- 
goliubov modes could be considered analogous to a recent 
proposal for topological phononic crystals, where sound 
waves propagate on top of a lattice of flowing air m- 
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Conclusion. —We proposed a generic scheme to create 
topological Bogoliubov excitations in systems described 
by a 2D nonlinear Schrodinger equation — a univer¬ 
sal equation governing systems ranging from exciton- 
polariton condensates to cold atom condensates and non¬ 
linear optical media. The key to our proposal is a back¬ 
ground mean field exhibiting vortex-antivortex pairs and 
a Kagome intensity pattern, which can typically be gen¬ 
erated by optical means. By virtue of weak interactions, 
Bogoliubov excitations propagating on top of this con¬ 
densate acquire non-trivial phases which break TR sym¬ 
metry and grant access to topological states. The re¬ 
sulting topological Bogoliubons manifest as robust chi¬ 
ral edge states, as we have explicitly demonstrated us¬ 
ing both tight-binding and full wave-expansion meth¬ 
ods. Our scheme does not require any external mag¬ 
netic field or net transfer of orbital angular momentum, 
which allows to generate chiral edge modes while keep¬ 
ing spin/polarization degrees of freedoms degenerate and 
simultaneously accessible. 
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